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ABSTRACT 


In  this  report,  ^  approximate  method  of  solution  is  given  to  determine  the 
location  of  the  rear  stagnation  point  in  the  viscous  base  flow  region  together  with 
the  size  of  the  viscous  wake  and  the  flow  characteristics  there.  It  is  useful 
as  a  means  of  obtaining  the  initial  conditions  for  the  near  wake  calculations.  The 
results  obtained  show  that  one  must  consider  the  finite  boundary  layer  at  the  base 
in  determining  the  location  of  the  rear  stagnation  point  and  the  viscous  wake  width 
there.  In  general,  the  viscous  wake  at  the  rear  stagnation  point  consists  of  a 
shell  of  hot  gas  surrounding  a  cool  core. 
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SYMBOLS 


X,  r  axial  and  radial  coordinates,  respectively,  non-dlmensionalized  with 

respect  to  6  =  „ 

D  the  density  non-dimensionalized  with  respect  to 

u,  V  the  axial  and  radial  velocity  components  non-dimensionalized  with 

respect  to  u 

o 

/  2 

p  the  static  pressure,  non-dimensionalized  with  respect  to  1/2  p  u 

o  o 

h  the  enthalpy,  non-dimensionalized  with  respect  to  h^ 


'B 
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B 


R 


eo 


R 


B 
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P 

r 


the  base  pressure 

the  free  stream  pressure 


the  base  wall  enthalpy,  non-dimensionalized  with  respect  to  h 


the  body  wall  temperature 

c  u  6 

OOP 

O 


the  base  radius 
2 


2  h 

o 


the  distance  along  the  body  surface  from  the  stagnation  point  to  the 
corner  of  the  base 


the  Prandtl  number 


6  the  viscous  wake  widtti  (  or  radius),  non-dimensionalized  with  respect 

to  6 

o 


iv 


B 


09 


n 


the  boundary  layer  thickness  at  the  corner  of  the  base 


the  coefficient  of  viscosity  non-dimensionalized  with  respect  to  U 
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the  expansion  angle 
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as  defined  in  equation  (15) 
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SUBSCRIPTS 


e  inviscid  condition  at  the  edge  of  the  viscous  wal<e 

o  condition  at  x  ==  o 

1  condition  at  the  rear  stagnation  point  x  =  xj^ 
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THE  APPROXIMATE  FLOW  CHARACTERISTICS  IN 


THE  BASE  REGION  OF  A  HYPERSONIC  AXI-SYMMETRIC  BODY 

INTRODUCTION 

A  theory  of  base  flow  for  h3T3ersonic  axi-symmetric  bodies,  blunt  or 
slender,  has  been  given  in  Reference  1.  It  was  developed  with  the  assumption 
that  the  gas  is  in  thermo-chemical  equilibrium  with  the  Lewis  number  for  the 
individual  chemical  species  equal  to  one.  The  integral  method  was  then  used 
to  reduce  the  governing  partial  differential  equations  into  ordinary  differential 
equations.  The  solutions  of  the  resulting  equations  and  the  appropriate  boundary 
conditions  yielded  the  detailed  flow  characteristics  in  the  base  region. 

In  this  report,  an  approximate  method  is  presented  by  which  the  location 
of  the  rear  stagnation  point,  the  size  of  the  viscous  wake,  and  the  flow  characteris¬ 
tics  there  can  be  determined  readily.  It  is  developed  based  on  the  analysis  given 
in  Reference  1  and  the  following  observations. 

It  is  noted  first  that  the  integi'al  continuity,  momentum  and  energy  equations 
(Equations  (1),  (2),  and  (3)  respectively)  can  be  integrated  further  when  the 
functional  form  of  ^  (x)  is  known.  From  available  shadowgraph  pictures  on  wake 
flow  in  the  base  region,  the  size  of  the  viscous  core  is  seen  to  decrease  almost 
conically  downstream  from  the  base  up  to  the  neck  region.  Therefore,  as  a 
first  approximation,  one  may  assume  that  6  (x)  is  a  linear  function  of  x.  Hence, 
the  three  equations  mentioned  previously  can  be  reduced  to  algebraic  equations. 

Furthermore,  with  the  occurrence  of  the  recirculating  flow  in  the  base 
region,  a  rear  stagnation  point  exists  which  is  imbedded  in  the  viscous  region 
and  along  the  axis.  Upstream  of  this  point  is  the  recirculating  flow  region.  The 
flow  in  the  mixing  layer  above  the  recirculating  flow  will  go  downstream  of  this 
point  and  form  the  neck  and  beginning  of  the  near  wake.  At  this  point,  the  velocity 
component  u  must  also  be  equal  to  zero  by  definition  since  v  is  already  equal  to 
zero  there.  This  condition  reduces  the  differential  equations  balancing  the 
momentum  and  energy  along  the  axis  (Equations  (5)  and  (6)  )  into  algebraic 


1 


relations.  In  this  manner,  a  whole  set  of  algebraic  equations  can  be  obtained, 
the  solutions  of  which  determine  the  location  of  the  rear  stagnation  point  and  the 
viscous  wake  size  there,  together  with  the  flow  characteristics  parameters. 

The  advantages  offered  by  the  present  method  are  manifold.  The  solution 
of  a  set  of  algebraic  equations  can  be  more  readily  obtained  than  the  original 
differential  equations.  When  the  external  inviscid  flow  conditions  are  uniform 
(the  so-called  "free-stream  line"  case),  a  table  can  be  generated  once  and  for 
all  to  give  the  solutions  in  terms  of  parameters  covering  a  wide  range  of  flight 
conditions  and  for  different  body  geometries.  Furthermore,  the  determination 
of  the  appropriate  boundary  for  the  inviscid  flow  field,  hence  the  base  pressure, 
requires  an  iteration  procedure  between  the  solutions  of  the  viscous  core  and  the 
inviscid  flow.  To  accomplish  this  would  be  a  lengthy  and  time-consuming  process 
if  the  theory  given  in  Reference  1  was  used  repeatedly.  This  is  especially  true 
for  a  slender  body  where  the  base  region  is  expected  to  be  long.  Therefore,  one 
can  use  the  method  and  results  given  in  this  report  and  arrive  at  a  reasonable 
boundary  shape  for  the  inviscid  flow  field  and,  hence,  the  base  pressure  quickly. 
Satisfactory  results  for  the  details  of  the  viscous  flow  field  in  the  base  region 
can  then  be  obtained  using  the  complete  theory  without  further  iteration. 

For  a  blunt  body,  the  base  region  is  generally  short  compared  with  the 
near  and  far  wake.  In  this  case,  the  detailed  knowledge  of  the  flow  characteristics 
in  the  base  region  may  not  be  important  with  respect  to  the  observables.  The 
present  method  then  gives  a  good  initial  flow  condition  for  the  near  and  far  wake 
calculations  with  the  effects  of  the  boundary  layer  and  the  subsequent  free  mixing 
being  accounted  for. 

For  a  chemical  flow  system  having  a  large  number  of  chemical  species,  the 
number  of  equations  to  be  treated  will  be  greatly  increased.  The  understanding 
of  such  a  complex  problem  is  limited  thus  far  mainly  because  of  the  mathematical 
difficulties  presented  in  solving  the  large  number  of  differential  equations.  It  is 
hoped  that  the  present  method  of  approach  can  be  extended  to  treat  such  a  problem 
for  the  base  flow.  If  the  governing  equations  again  can  be  reduced  to  algebraic 
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equations  in  the  same  manner  as  in  the  present  case,  it  is  possible  that  solutions 
can  be  obtained  in  the  neck  region.  The  results  should  serve  as  a  good  guide  in 
developing  a  more  complete  and  detailed  theory. 
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1.  THE  GOVERNING  EQUATIONS 

For  convenience,  the  governing  equations  given  in  Reference  1  which  are 
pertinent  to  the  present  analysis  are  summarized  here. 

The  continuity,  momentum  and  energy  equations  after  their  integrations 
across  the  wake  are 
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(cont) 


Two  more  ordinary  differential  equations  are  obtained  by  satisfying  the 
momentum  and  energy  balances  along  the  axis  r  =  0.  They  are: 
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77  =  o.  The  outer  boundary  conditions 
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give  two  additional  algebraic  equations. 

The  profiles  are  assumed  to  be 

2 

~p~^  "  ^7)  +  (x)  (1  -  +  a^  (x)  (1  -  Tj^)  +  a^  (x)  (1  -  r\^)^ 

(9) 

-^  =  G  (x,  rj)  ^  b^  (X)  (1  -  if)  ,  b^  (X)  (1  -  rf)  +  bg  (X)  (1  - 

(10) 

where  the  a's  and  the  b's  are  the  unknown  quantities.  They  satisfy  the  boundary 
conditions  that 


(1)  at  ?7  =  o, 

(2)  at  Tj=  1, 

(3)  at  X  =  o 


u  =  u  V  =  V 
e  e 


pu  =  F  (o,  r) 


h  =  G  (o,  r) 


=  o 


and  p  =  p 


(11) 


(12) 


(13a) 

(13b) 


where  F  (o,  r)  and  G  (o,  r)  are  the  given  and  known  profiles  at  the  base  (i.e. , 
^  ~  o)*  The  initial  profiles  F  (x,  tj)  and  G  (x,  t])  have  been  assumed  as 


F(x,  rj) 
G(X,  Tj) 

F(x,  Tj) 

G(x,  tj) 


(14a) 


(14b) 


6 


where 


= 


6^ 

6(x) 


(See  Figure  1) 


In  here,  h^  is  the  assumed  constant  nondimensionalized  wall  enthalpy  at 
the  base.  The  fimctions  f  and  g  have  the  same  functional  forms  as  the  distribu¬ 
tions  of  the  pu  and  h  across  the  boundary  layer  at  the  corner  of  the  base  of  the 
body,  respectively.  The  function  77^  was  assumed  as 


^d  = 


6(x) 


=  1- 


=  1- 


6(x) 

6(x) 


1+2  _55e_  _x_ 

^e"e  "b 


1  +  0, 


i) 


for  laminar  flow 

(15a) 


for  turbulent  flow 
(15b) 


where  6  is  the  known  boundary  layer  thickness  at  the  base 

s  is  the  distance  along  the  body  surface  from  the  stagnation  point  to 
the  corner  of  the  base. 

The  parameter  a  (or  )  characterizes  the  relative  growth  of  a  laminar 
(or  turbulent)  free  mixing  layer  with  respect  to  that  of  the  laminar  (or  turbulent) 
boundary  layer  (Reference  1). 

Equations  (1)  to  (3)  and  (5)  to  (8)  form  the  basic  system  for  the  seven 

unknowns  6,  a  and  b  . 

s  s 


2.  THE  APPROXIMATE  METHOD  OF  SOLUTION 
(a)  The  General  Case 

It  is  noted  that  6(x)  may  be  assumed,  as  a  good  approximation,  to  vary 
linearly  with  x,  or 

1-61 

6=1-  -  X  (16) 

^1 
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where  6^  is  the  unknown  nondimensionalized  wake  radius  at  the  unknown  location 
X  =  of  the  rear  stagnation  point.  With  this  assumed  functional  form  (16)  for 
6(x),  Equations  (1),  (2)  and  (3)  can  be  integrated  from  x  =  0  to  x  =  x^^  to  yield 
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The  right-hand  sides  of  Equations  (17),  (18)  and  (19)  are  all  integrable  as  the 
functions  p  ,  u  ,  v  ,  h  ,  ^  i  and  Q  are  all  given  functions  and  6(x)  is 

assumed  as  in  (16).  One  exception  is  the  function  4-  which  is  the  dissipation 
existed  in  the  flow.  In  the  base  region,  the  contribution  to  the  dissipation  due  to 
the  recirculating  flow  is  exijected  to  be  small.  Hence,  as  an  approximation,  one 
may  consider  only  the  dissipation  from  the  boundary  layer  flow  and  the  subsequent 
free  mixing  in  the  calculation  of  the  function  4.  In  other  words,  4  will  be  cal¬ 
culated  by  considering  only  the  terms  F(x,  tj)  and  G(x,  rj)  in  the  assumed  pro¬ 
files.  Therefore,  4  is  also  a  known  function  of  x. 

Equations  (17),  (18)  and  (19)  are  three  algebraic  equations  relating  the 
known  inviscid  flow  conditions  to  the  unknowns  6^,  x^,  and  other  unknowns  in  the 
assumed  profiles.  Additional  algebraic  equations  are  now  obtained  by  considering 
the  conditions  at  the  rear  stagnation  point. 


By  definition,  at  the  rear  stagnation  point  x  =  x^  and  r  =  o,  the  velocity 
component  us  must  be  equal  to  zero  since  v  =  o  at  r  =  o  in  general.  Therefore,  from 
(9)  one  has 
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The  last  two  equations  can  be  reduced  to,  from  the  assumed  profiles  (9)  and  (10), 


/  n  J  V  2  (u  u  ),  r 

(•  T  ■&),  *  J  I' 


1  R  6 
®o  1 


(-2aii  -  4a2j  -  Oa^^) 
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(22) 


(23) 


where  a  ,  a  ,  a  ,  b  ,  b  ,  and  b  are  their  respective  values  at  x  =  x  . 
XX  ZX  oX  XX  ^X  oX  X 

Furthermore,  a^^  and  b^^^  should  satisfy  the  following  relations  through  (7)  and 
(8); 


(24) 
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(25) 


where 


Equations  (17)  to  (19)  and  (21)  to  (25)  are  the  eight  algebraic  equations  to 

be  solved  simultaneously  for  the  eight  unknowns,  x,,  6,»  a  and  the  b  .  In  the 

11s  s 

general  case,  they  can  be  reduced  into  three  algebraic  equations  for  x^,  and 
a  .  When  the  external  velocity  and  enthalpy  gradients  are  negligible,  i.e. , 

a  X 

F^',  ,  uu  and  Q  are  approximately  equal  to  zero,  a^^  and  bj^^  are  equal  to  zero 

according  to  equations  (24)  and  (25).  The  remaining  equations  can  then  be  re¬ 
duced  to  two  equations  for  x^^  and  6^  in  such  a  case. 
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(b)  The  "Free-Stream  Line"  Case 


For  the  particular  case  when  the  external  inviscid  flow  conditions  are 
uniform,  the  so-called  "free-stream  line"  case,  the  set  of  equations  given  in 
the  previous  section  offers  a  convenient  way  for  solutions.  In  this  case.  Equa¬ 
tions  (17)  to  (19)  and  (21)  to  (23)  can  be  rewritten  and  reduced  into 
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(26) 


(27) 


(28) 


(29) 


=*’21  ^  =’’31  '  "  <='> 

From  (29)  and  (30),  one  has  a  =  a  =0  and  from  (31)  b  =  (-2/3)b  . 

uX  ^X 

Therefore,  there  are  only  three  unknowns,  6,  x  ,  and  b  ^  left  to  be  solved  from 

Equations  (26)  to  (28).  It  can  be  shown  that,  in  general,  ©  and  ©i,  are  linear 

1 

fimctions  of  b  .  Therefore,  b  can  be  eliminated  simply  by  combining  Equa- 
tions  (27)  and  (28).  The  resulting  equation  together  with  Equation  (26)  deter¬ 
mines  and  A  convenient  method  of  solution  for  these  equations  can  be 
devised  such  that  tables  can  be  generated  once  and  for  all  to  give  solutions  in 
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terms  of  parameters  covering  a  wide  range  of  flight  conditions  and  for  different 
body  geometries.  This  will  be  explained  in  detail  in  the  following  section. 

3.  NUMERICAL  PROCEDURES 

For  illustrative  purpose,  let  us  consider  the  "free  stream-line"  case  and 
assume  that 
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for  laminar  flow 
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for  turbulent  flow 
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Therefore,  Equations  (27)  and  (28)  can  be  combined  to  eliminate  b  .  This 
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V, ‘®<>  ■  ®11  «/' -  821 

21  o  11 


=  <8h^,  -  821 


'•i 


2m  CD 

"  B  821  ’‘l 


(35) 


In  principle,  then,  6^,  and  x^^  can  be  obtained  by  solving  Equations  (26) 
and  (35).  However,  is  a  function  of  both  the  unknowns  x^^  and  for  a  given 
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condition  j^Equation  (34)J .  Therefore,  to  obtain  the  solution  of  the  simultaneous 
equations  for  6^  and  directly  will  require  a  somewhat  lengthy  and  laborious 
procedure.  But  a  convenient  way  for  solutions  can  be  obtained  in  the  following 
manner. 

Let  us  rewrite  (26)  to  give 


ilii  ^  ^  ,.2. 

X,  oil  2  I  1  ^o  1  "1 


(36) 


and 


2  X 

2(^o'^i  ^i"  “i - ^ 


Vo(l  +  6j) 


Putting  these  into  Equation  (35),  one  has 


(1  +  6,)  6,“  Tejj  0 

^  A  X  J.  X 


'”21  0  -■ 


_4m 

V 


ncD  r  2  1  1  1 


o  e 


(37) 


This  equation  gives  6^^  as  a  function  of  h^,  m,  v^,  and  A^^  where  Aj^  is 

a  function  of  6  (=  A  ),  a{or  a  )>  s„,  x  and  6  .  For  a  given  body  geometry  at 
Ho  f  H  1  1 

given  flight  conditions,  m,  Rp  ,  6t,(=  A  ),  s  are  known  quantities  and  for  an 
assumed  expansion  angle  around  the  corner  of  the  base  (or  assumed  base  pres¬ 
sure),  h  ,  V  and  ft  (or  rr  )  can  be  determined.  To  solve  for  6,  and  x  one  may 
B  o  f  11 

assume  a  value  of  A  which  varies  from  A^  1°  o"®  solve  for  6^^, 

from  (37)  for  a  given  set  of  values  of  h_,  m,  Rp  ,  v  ,  A  »  a{orCTJ 

B  o  o  o  i 

s  .  X  is  then  determined  from  (36).  For  this  set  of  values  of  x  and  6  , 

HI  1  -L 

A,  can  be  calculated  from  (34)  for  the  same  given  values  of  6„,  a  and  s  .  This 
1  H  H 

calculated  value  of  Aj^  should  check  with  the  assumed  Aj^  value  when  6^^  and  x^  are 
the  correct  solutions  to  both  Equations  (35)  and  (26). 
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However,  this  trial  and  error  procedure  is  not  really  necessary.  Instead, 

one  may  solve  for  6^  from  (37)  and  then  from  (36)  assuming  a  series  of  values 

of  X,  from  X  to  one  for  a  given  set  of  values  h„,  m,  Rp  ,  v  ,  X  and  s„.  rr 
1  o  B  o  0  o  B 

and  a  are  then  calculated  from  |^see  (34)J  . 


a 


(38a) 


a 

( 


^6 

v 


(38b) 


In  this  manner,  tables  can  be  generated  beforehand  to  give  the  values  of  6^,  x^, 

CT  and  cr ,  as  functions  of  X,  for  a  given  set  of  values  of  h.,,  v  ,  X  ,  s„,  m  and 
f  1  c  o  o  B 

Re^.  For  a  particular  case  of  interest,  one  merely  has  to  go  to  these  tables  and 

find  out  the  values  of  6,  and  x,  for  a  particular  set  of  values  of  1l  ,  v  ,  X  i  s  , 

11  B  o  o  B 

m,  Rg^  and  a  (or  )•  The  corresponding  value  of  Xj^  which  gives  the  ratio  of 
the  thickness  of  the  mixing  layer  to  the  viscous  wake  width  at  the  rear  stagnation 
point  can  also  be  obtained. 

To  summarize,  the  numerical  procedure  consists  of:  (1)  construction  of 
tables  (or  graphs)  giving  values  of  a.  (r  6^.  x^,  and  b^^  (from  Equations  (38a), 
(38b),  (37),  (36),  and  (27)  respectively)  as  functions  of  Xj^  for  sets  of  assigned 
values  of  s  ,  >  (=  g  ),  m,  Rp  ,  v  ,  and  h  .  (2)  For  a  particular  case  of  interest, 
s  is  known  from  the  given  body  geometry.  X  (=  6-r,).  ni,  and  Rg  can  be  obtained 
from  the  known  or  calculated  inviscid  shock  layer  flow  and  the  boundary  layer  flow 
for  a  given  flight  condition.  (3)  Assume  a  base  pressure  value  hence  the  value  of 
the  expansion  angle  Q  and  determine  the  values  of  v  =  -  tan  6 ,  h^  and  a  (or  a  )• 
(4)  The  solutions  for  6  ,  x  ,  and  b  are  then  obtained  from  the  tables  for  this 

X  X  ^  J. 

set  of  values  of  s  ,  X  .  m,  Rg  ,  v  ,  h„  and  q-  (or  cr  ). 

B  O  ^O  O  B  C 

It  should  be  remarked  that  Xj  varies  between  X^  and  one  only  since  77 ^  = 
6di/6i  can  not  be  negative  [see  Equation  (34)J  .  Also  the  value  of  a  (or  (t^)  is 
greater  than  one  in  actual  case  because  the  growth  of  a  free  mixing  layer  is 
larger  than  that  of  a  boundary  layer  under  the  same  condition. 


15 


4.  RESULTS  AND  DISCUSSIONS 

Results  given  in  this  section  are  for  the  "free-stream  line"  case  with 
F(x,  Tj)  and  G(x,  tj)  given  in  (32a)  and  (32b). 

Tables  1(A)  and  (B)  give  part  of  the  results  obtained  in  the  manner  described 

in  step  (1)  of  the  preceding  section.  In  these  calculations,  the  term  on  the  right 

hand  side  of  (37)  has  been  found  to  be  small  and  hence  neglected.  These  two 

tables  are  obtained  for  s„  =  4,  v  =  -0.06993. ,  h_  =  2  and  1/2  with  X  =0.1 

Bo  B  o 

for  table  1(A)  and  X  =0.001  for  table  1(B).  The  range  of  a  values  covered  is 
o 

approximately  from  3.2  to  1. 

Figures  2a  and  2b  are  plots  of  the  results  given  in  Tables  lA  and  IB, 

respectively.  They  show  that  6^^  decreases  and  x^  increases  with  cr  except  that 

increases  slightly  beyond  (j  =  2.8  in  the  case  of  =  0. 1.  Figure  2a  shows 

further  that  6,  increases  and  x,  decreases  with  h^  in  the  case  of  X  =0.1  but 
1  1  B  o 

they  are  practically  unchanged  with  h  for  >  =  0.001.  However,  it  should  be 
remembered  that  v^,  h^  and  (j  are  not  independent  of  each  other  but  rather  they 
are  related  through  the  base  pressure  (or  the  expansion  angle)  assumed.  There¬ 
fore,  one  should  not  draw  any  definite  conclusion  from  results  such  as  in  Table 
lA  (or  Figure  2A)  in  regard  to  the  variations  of  and  x^^  with  respect  to  a  or 
hg  Independently  for  a  particular  given  case  since  neither  a  nor  hg  is  a  single 
independent  variable  parameter.  This  will  be  further  explained  later  in  this 
section. 


However,  one  can  draw  a  conclusion  as  to  the  effect  of  the  finite  boundary 

layer  at  the  base  from  these  results.  Figure  3  is  a  plot  of  6^^  and  x^^  vs  cr  for 

6  =  0. 1  and  0.001  with  the  same  values  of  s  ,  v  and  11.  Since  for  a  given 

B  B  o  H 

body  geometry,  flight  conditions  and  an  assumed  expansion  angle,  the  values  of 
®B’  '^o’  ^  ^  known  and  fixed,  therefore,  this  figure  com¬ 

pares  the  results  of  6  and  x  from  the  "actual"  case  (6  =0.1  case)  with  the  case 

-I-  B 

where  6g  is  assumed  negligible  (6g  =  0.001  case).  It  is  seen  that  if  fig  is  neglected, 
is  much  smaller  and  x^^  much  larger  than  the  actual  case.  The  discrepancy 
increases  with  increasing  (j  and  varies  from  15%  to  50%  for  and  30%  to  50% 
for  Xj^.  This  can  be  explained  as  follows.  For  an  initial  boundary  layer  of  small 
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but  finite  thickness  at  the  base,  the  volume  of  gas  =  1  -  (1  -  6g)^J  contained 
in  the  boimdary  layer  there  is  not  negligibly  small.  In  the  axi-sjnn metric  case, 
the  total  volume  of  this  gas  surrounding  the  base  converges  to  the  axis  to  form 
the  core  at  the  rear  stagnation  point.  The  radius  of  this  core  must  be  sufficiently 
large  so  that  this  finite  volume  of  gas  will  be  contained  there.  When  the  initial 
boimdary  layer  at  the  base  is  neglected,  under  the  same  conditions,  the  core  size 
will  be  naturally  much  smaller  since  there  will  be  very  little  volume  of  gas  to 
be  contained  in  the  core  in  this  case.  At  the  same  time,  will  be  larger  for 
the  e:q>ansion  angle  (or  base  pressure).  With  a  larger  value  of  (7,  more  volume 
of  gas  will  be  drawn  from  outside  into  the  mixing  layer  and  hence,  the  difference 
in  and  x^^  will  become  larger  between  the  two  cases.  These  results,  therefore, 
emphasize  the  importance  of  the  inclusion  of  the  finite  boimdary  layer  in  treating 
the  base  flow  for  an  axi-symmetric  body. 


Table  n  gives  results  for  the  specific  case  of  a  20°  sphere  cone  at  = 

19.52  and  altitude  =  150,000  feet.  They  are  obtained  according  to  the  procedure 

described  in  the  preceding  section.  In  these  calculations,  two  constant  body  wall 

temperature  values  Tw_  (2  and  1/2  corresponding  to  hot  and  cool  wall)  have  been 

assumed  and  the  6  values  were  estimated.  In  each  case,  five  expansion  angles 
B 

and  their  respective  base  pressure  values  were  assumed  and  the  corresponding 
sets  of  values  of  v^,  h^  and  (j  were  determined.  The  solution  for  6^,  x^^,  b^^^ 
and  the  value  of  A  ^  were  then  obtained  from  the  table. 

These  results  show  that  5^  decreases  and  x^  increases  with  increasing 
Twg  ^  the  same  base  pressure  or  with  increasing  base  pressure  for  the 
same  T^^  or  h^.  However,  it  would  be  clearer  and  more  meaningful  if  one  uses 
the  parameters  x^^  and  0^^  for  comparison  instead  of  x^  and  6  where  0^^  denotes 
the  angular  position  of  the  edge  of  the  viscous  wake  with  respect  to  the  free 
stream  (see  Table  II  and  Figures  4  and  5).  These  show  that  0^  decreases  and 
Xj^  increases  with  increasing  base  pressure  for  the  same  T^^^  or 
smaller  0^^  and  larger  x^^  will  give  a  larger  base  region,  therefore,  the  larger  the 
assumed  base  pressure,  the  larger  is  the  base  region.  This  is  expected  phys¬ 
ically  since  a  higher  base  pressure  can  sustain  the  mixing  process  longer  and  a 


Since  a 
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bigger  volume  of  gas.  Furthermore,  is  practically  unchanged  but  increases 
with  h  for  a  given  base  pressure.  This  is  again  because  a  hot  base  wali  can 
sustain  a  longer  mixing  process.  Furthermore,  for  a  hotter  base  wall  but  with 
the  same  base  pressure,  the  resulting  mixing  layer  at  the  rear  stagnation  point 
is  also  thicker  as  it  is  evident  from  the  larger  value  of  X  ^  and  smaller  value  of 

This  has  the  same  trend  as  in  the  case  of  the  boundary  layer.  However, 
it  is  noted  that  the  base  pressure  and  the  base  wall  enthalpy  are  not  two  inde¬ 
pendent  parameters.  For  a  given  known  base  wall  enthalpy,  there  is  only  one 
corresponding  base  pressure.  This  will  be  a  subject  for  future  investigation. 

The  other  interesting  results  are  those  of  b^^  and  6(j^.  The  value  of  b^^^ 
is  practically  zero  for  all  the  cases  considered.  This  means  that  the  radial 
enthalpy  distribution  at  the  rear  stagnation  point  in  the  viscous  wake  is  the  same 
as  that  at  the  base  but  in  different  scale  (see  Figure  5).  It  is  further  noted  that 
6dj^  is  generally  nonzero  at  x^^.  In  the  case  of  a  cool  base,  this  means  that  a 
cool  core  exists  there  with  a  shell  of  hot  gas  surrounding  it.  This  shell  of  hot 
gas  is  formed  mainly  by  the  gas  from  the  boundary  layer.  For  a  hot  base,  the 
solution  shows  that  the  core  is  hot.  However,  in  the  actual  case  the  temperature 
and  enthalpy  at  the  wall  of  the  base  probably  will  vary  rather  than  be  uniform  as 
assumed  here  and  decrease  from  the  value  at  the  corner  to  some  lower  value  at 
the  center.  If  this  is  the  case,  then  even  for  an  originally  hot  body  wall,  that 
part  of  the  base  near  the  axis  is  cool.  Therefore,  at  the  rear  stagnation  point, 
the  core  will  still  be  cool  but  the  shell  surrounding  it  will  be  hot  and  compara¬ 
tively  thicker  than  in  the  case  of  a  cool  wall.  These  phenomena  are  probably 
associated  with  the  growth  of  the  mixing  layer.  Over  the  extent  of  the  base 
region,  up  to  the  rear  stagnation  point,  there  would  not  be  sufficient  time  or  dis¬ 
tance  for  the  mixing  layer  to  grow  to  its  fullest  before  its  lower  boundary  reaches 
the  axis.  The  heat  content  is  being  contained  within  the  mixing  layer  without  too 
much  heat  exchange  between  the  mixing  layer  and  the  recirculating  flow.  There¬ 
fore,  the  enthalpy  level  between  the  axis  and  the  lower  boundary  of  the  mixing 
layer  would  be  the  same  as  that  of  the  base  wall  with  the  enthalpy  distribution 
in  the  mixing  layer  remaining  the  same  as  in  the  initial  boundary  layer.  How¬ 
ever,  future  experimental  measurements  and  flight  tests  for  the  base  region  will 
shed  additional  light  on  this  problem. 
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It  should  be  noted  that  the  location  of  the  rear  stagnation  point  does  not  neces¬ 
sarily  coincide  with  the  location  of  neck  which  is  the  narrowest  region  of  the 
wake.  However,  the  results  obtained  for  the  rear  stagnation  point  are  necessary 
for  the  investigation  of  the  flow  in  the  near  wake  region. 

5.  CONCLUSIONS 

It  can  be  concluded  from  the  results  in  the  present  investigation  that  one 
must  consider  the  finite  boundary  layer  at  the  base  in  treating  the  base  flow.  By 
neglecting  the  boundary  layer  one  will  overestimate  and  underestimate  6^^  by 
a  large  percent. 

For  an  uniformly  cool  base,  it  was  found  that  a  cool  core  exists  at  the 
rear  stagnation  point  with  a  shell  of  hot  gas  surrounding  it.  This  shell  of  hot  gas 
is  formed  mainly  by  the  gas  from  the  boundary  layer.  However,  a  hot  core  could 
be  formed  there  if  the  base  were  uniformly  hotter  than  the  free  stream.  The 
latter  case  may  not  be  realistic  since  the  base  wall  temperature  will  probably 
decrease  from  the  corner  of  the  base  to  the  axis.  Therefore,  it  is  reasonable 
to  say  that,  in  general,  a  shell  of  hot  gas  surrounding  a  cool  core  exists  at  the 
rear  stagnation  point. 
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TABLE  I 


W  ^  =  'b  =  “-I’  " 

V  =  -  Tan  0  -  -0.  06993 


(a)  hg  =  2 


a 

a 

e 

6 

1 

X 

1 

’^21 

S 

3. 1694 

1.78 

0. 4629 

8. 1324 

-0.0848 

■ 

2.8137 

1.545 

0.4584 

8.0949 

-0.0767 

2. 5080 

1.34 

0.4603 

7.9881 

-0.0625 

2. 2381 

1.17 

0.4688 

7.8029 

-0.0451 

0.7 

1.9926 

1.01 

0.  4849 

7.5215 

-0.0286 

0.6 

1.7621 

0.872 

0.  5107 

7.1108 

-0.0161 

0.  5 

1.5385 

0.738 

0.  5505 

6. 5050 

-0.0079 

0.4 

1.3132 

0.604 

0.  6144 

5.  5606 

-0.0032 

0.3 

1.0741 

0.468 

0.7288 

3.8988 

-0.0009 

0.2 

(b)  hg  =  1/2 


a 

CT 

c 

6 

1 

X 

1 

•^21 

s 

2.8542 

0.4266 

8.4464 

-0.0164 

m 

2. 5906 

0. 4304 

8.3458 

-0.0160 

2.3493 

0.4384 

8. 1906 

-0.0141 

■H 

2. 1246 

1.085 

0.4516 

7. 9683 

-0.0107 

0.7 

1. 9113 

0.956 

0,4712 

7.6579 

-0.0070 

0.  6 

1. 7045 

0.83 

0. 4999 

7,2230 

-0.0039 

0.  5 

1. 4990 

0.707 

0. 5422 

6. 5948 

-0.0019 

0.4 

1.2877 

0.  587 

0.6085 

5. 6274 

-0.0007 

0.3 

1. 0594 

0.457 

0.7255 

3.9390 

-0.0002 

0.  2 
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TABLE  I  (Continued) 


(B)  X  =  6  =  0.001,  S_  =  4 

O  t>  -D 


V  =  -  Tan  0  =  -  0,  06993 
o 


(a)  lig  "  2 


a 

a 

e 

6 

1 

X 

1 

^1 

X 

1 

3. 6284 

0.1599 

12. 1182 

0 

3. 1660 

0.1841 

11.  7690 

2. 6217 

0. 2248 

11. 1814 

0.02 

1.9170 

0.3170 

9.8511 

0.01 

1. 7377 

0.788 

0.3542 

9.3143 

0.008 

1. 5340 

0.68 

0.4088 

8. 5271 

0.006 

1.4192 

0.  62 

0. 4477 

7. 9661 

0.005 

1.2919 

0.  555 

0. 5004 

7.  2059 

0.004 

1. 1469 

0.483 

0. 5777 

6.  0916 

0.003 

0.9739 

0.39 

0.7073 

4. 2218 

0 

0.002 

(b)  hg  =  1/2 


a 

a 

c 

S 

X 

1 

*^21 

3.6221 

0.1596 

12.1218 

0 

3. 1617 

0.1838 

11. 7720 

2. 6193 

0. 2246 

11. 1839 

0.02 

1.9161 

0.884 

0.3169 

9. 8528 

0.01 

1.7370 

0.788 

0.3541 

9.3157 

0.008 

1.  5336 

0.  68 

0.4087 

8. 5283 

0.006 

1.4188 

0.62 

0. 4476 

7. 9671 

0.005 

1.2916 

0.  555 

0. 5004 

7.  2067 

1 

0.004 

1. 1467 

0.483 

0. 5776 

6. 0922 

0.003 

0. 9738 

0.39 

0.  7073 

4. 2222 

0 

0.002 
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FIGURE  |:  FLOW  IN  THE  BASE  REGION 
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FIGURE  5:  ENTHALPY  DISTRIBUTIONS  IN  THE  BASE  REGION 
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In  this  report,  an  approximate  method  of  solution  is 
given  to  determine  the  location  of  the  rear  stagnation 
point  in  the  viscous  base  flow  region  together  with  the 
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obtained  show  that  one  must  consider  the  finite  boundary 
layer  at  the  base  in  determining  the  location  of  the  rear 
stagnation  point  and  the  viscous  wake  width  there.  In 
general,  the  viscous  wake  at  the  rear  stagnation  point 
consists  of  a  shell  of  hot  gas  surrounding  a  cool  core. 
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